Role of spinon and spinon singlet pair excitations on phase transitions in d — wave 

superconductors 
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We examine the roles of massless Dirac spinon and spin singlet pair excitations on the phase 
transition in d — wave superconductors. Although the massless spinon excitations in the presence 
of the spin singlet pair excitations do not alter the nature of the phase transition at T = 0, that 
is, the XY universality class, they are seen to induce an additional attractive interaction potential 
between vortices, further stabilizing vortex-antivortex pairs at low temperature for lightly doped 
high T c samples. 
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Recent thermal hall conductivity measurements 1 sug- 
gest the existence of vortices in the pseudogap (PG) 
phase. This implies that preformed pairs are present 
in the PG phase 2 . In the PG phase, vortex-antivortex 
pairs remain broken to cause a state of globally incoher- 
ent but locally coherent Cooper pairs. Vortex induced 
phase transitions in underdoped region have been an is- 
sue of great interest 3-8 . In this paper, by applying a dual- 
ity transformation we extend the U(l) gauge Lagrangian 9 
obtained from the slave-boson mean field theory 10,11 in 
order to examine how vortex induced phase transitions in 
d — wave superconductors at low energy are affected by 
the presence of both the massless spinon and spinon sin- 
glet pair excitations. In other study 7,8 , the flavor number 
of the massless Dirac fermions without the spinon singlet 
pair excitations is shown to alter the nature of the phase 
transition. According to this study 7 , as the flavor num- 
ber increases, the type II superconductivity is preferred 
showing the second order phase transition which deviates 
from the XY universality class. In the case of small flavor 
number leading to the type I superconductivity, it be- 
comes the first order transition owing to the strong fluc- 
tuations of the massless gauge field. Our present study 
differs from other previous studies 7,8 in that in our case 
the U(l) effective gauge field of interest becomes massive 
as a result of the spinon singlet pair excitations. Particle- 
hole excitations of the massless Dirac spinons lead to the 
renormalized kinetic energy of the U(l) gauge field 4 ' 12 ' 13 
(the second term in Eq. [3]). However, we find that the 
XY universality class is not altered despite the presence 
of the massless Dirac fermions (spinons). Thus it is ir- 
respective of the flavor number as long as the effective 
gauge particle (Eq. [5] and Eq. [6]) remains sufficiently 
massive. The U(l) Berry gauge field can emerge from 
coupling between the massless Dirac fermions and the 
vortices 4-6,8 , which may affect the XY universality class. 
However, the Berry gauge field becomes massive owing 
to the charge fluctuations 5 and as a consequence the ef- 
fective dual Lagrangian stays robust to maintain the XY 
universality class. It is shown from the present study that 
the interaction potentials between vortices are modified 



to bring an additional attractive interaction (Eq. [7] and 
Eq. [11]) as a consequence of the massive gauge field. 

Our primary focus is to examine how low energy exci- 
tations affect the vortex induced phase transition. Here 
low energy excitations refer to the phase fluctuations of 
both the spinon singlet pair and the single holon order 
parameter and the massless Dirac spinon excitations near 
the d—wave nodes of the spinon singlet pair. Gauge field 
fluctuations are introduced to allow the presence of inter- 
nal flux responsible for energy lowering. Thus considering 



proper phase fluctuations (involved with 



'sp 



and 



4>b = e * ) for the spinon pairing order field and the single 
holon order field respectively, we rewrite the low energy 
effective Lagrangian of Lee 9 in compact form, 
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is the renormalized 



where Kb iM = (l/ub,Kb,Kb) with 1/ub (~ 1/i), the 
compressibility and Kb (~ 2txoS), the phase stiffness 
of the single holon field and K sPifl = (l/u sp , K spi K sp ) 
with l/u S p (~ 1/t/), the compressibility and K sp (~ 
JAq), the phase stiffness of the spinon pair order 

field. ?/w = m ii n A 

Nambu spinor associated with the d — wave gap nodes 
n. vp (~ Jxo) and V A (~ JAq) are the fermi 
and gap velocities of the Dirac spinons respectively. 
J i> = \^ n ^\ia TZ ^o,iv F vb\ a vb lai iv F vb\ a vb 2a ) is the 
three current of the spinon quasiparticles and p sp , the 
average density of spinon pairs. 

By introducing the unitary gauge transformation a M = 
2a p + d^9 sp , we rewrite Eq. [1] 
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where 9 p = 9 sp — 29b- 9 p is the phase of the Cooper 
pair order parameter A Co oper(k) = |A^ ooper (/c)|</>p 
where A C oo P er(k) =< c fcT c_ fei >=< b*f k \b*f- ki >=< 
./• feT ./L fci X b* > 2 = |A (fc)|| < b* > \ 2 <l> sp( f>f and thus 
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fective gauge field associated with the phase fluctuations 
of spinon singlet pair order parameter and the original 
internal U(l) gauge field a M . The mass of the effective 
gauge field a M is defined by the phase stiffness K sp of the 
spinon pairing order parameter associated with the PG 
phase of the doped Mott insulator. In the above equa- 
tion, the fluctuating fields of 9 P , a M and tp n , are U(l) 
gauge invariant, thus satisfying the Elitzur's theorem 14 . 

Integrating over the Nambu spinor fields and expand- 
ing the resulting logarithmic term upto second order in 
we get an effective U(l) Lagrangian involved 



= 4,12,13 



a 



with the massive gauge field a M , 
Z = J Dd p D~a^e-V dT f dx2c , 
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where Kb lfl = K^^/A and m 2 „ = K spp . N is the number 
of flavors (i.e., the number of nodal points) of the Dirac 
fermions. The kinetic energy term (the second term) 
of the effective gauge field a p arises as a result of the 
massless excitations of the spinon quasiparticles (Dirac 
fermions) 4,12,13 . The Berry phase contribution ip sp d T 9 p 
is related to the Cooper pair boson density 15 ' 16 . 

If we ignore the spinon quasiparticles and thus con- 
sider only the N = limit (which corresponds to the 
isotropic s — wave superconductivity) , it is obvious that 
the kinetic energy term of the gauge field disappears. In- 
tegrating over the effective gauge field a M in Eq. [3], we 
obtain 



K 



d^Op + 2A P \ 2 + ip sp (d T p + 2A ) 
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where K v 



K b ^K 3 



This expression shows that 



the phase stiffness K Pi{1 of the Cooper pair field is in a 
reduced "mass" form and the usual logarithmic type of 
interaction between vortices arises. 

The duality transformation 3,5,15-17 of Eq. [3] with the 
introduction of vortex mass and self interaction terms 
leads to an effective Lagrangian for the vortex field, 
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where m 2 v ~ Kb — Kb c ~ 6 — S c with S, the hole doping 
concentration and S c , the critical hole doping concentra- 
tion and a — -S^- — Ubpsv 4> P v represents the Cooper 

pair vortex field and c M , the dual gauge field to mediate 
interactions between vortices. Noting that p, is analo- 
gous to an applied "magnetic field" H z and (d x c) T to 
-Bz 15,16 , the sixth term in Eq. [5], —p(d x c) T 15,16 which 
results from the Berry phase term ip sp d T 9 p is analogous 
to interaction energy —H Z B Z associated with the vortex 
field. We note that in case of 5 < S c vortex condensation 
occurs. 

Integrating over the effective gauge field o M , we get 



Z = Z« J D^vDc^e- J d3xC > 
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-idx c-2A-p{dx c) T , (6) 



The fifth term represents an additional kinetic energy of 
the dual gauge field resulting from coupling of the mass- 
less Dirac spinon field to the Cooper pair field via the 
massive effective gauge field o M . Not considering the ex- 
ternal magnetic field p, despite the contribution of the 
massless Dirac fermions (as shown in the fifth term) the 
nature of the XY universality class will not be affected 
as long as the mass of the gauge field, that is, the phase 
stiffness K sp of the spinon singlet pair order parameter 
is substantially large; the lower temperature, the larger 
K sp in the PG phase. 

By considering a static case of the vortex field we cal- 
culate the dual gauge field propagator to obtain the in- 
teraction potential between vortices. The dual field prop- 
agator involved with the fourth and fifth terms in Eq. [6] 
is obtained, 
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where K p = M = >/M K 'P + K b/*) and 

zj = Kb +4 K ■ i,j = 1,2 denotes the space component 
index (x,y). In the real space, the above equations lead 
to the total interaction potential between a vortex and 
an antivortex, 



2 



V(x-x') = 2nK p ln\x- x'\ 

-~Tr 2 K b zj(StruveH (M 2 \x - x'\) - Y (M 2 \x - x'\)) 
— 2irK p ln\x - x'\ + n 2 K b z,jln\x - x'\, (8) 
for M 2 \x-x'\ « 1 



~ 2nK p ln\x -x'\- ^fi' J 
for M 2 b-a;'| >> 1 



x — x' 



(9) 



where Yo(x) is the zeroth order Bessel function of the 
second kind and StruveHo(x) is the zeroth order Struve 
function 18 . The additional interaction between vortices 
shows a power law decay at large separation distances be- 
tween vortices while it is logarithmic in nature at short 
separations. In Fig. 1, we plot this interaction potential 
as a function of distance in the underdoped region 19 . The 
correction term yields a large modification to the loga- 
rithmic interaction energy particularly at short distances, 
showing additional attractive interactions between vor- 
tices. It becomes more effective at lower temperature 
and causes stronger vortex-antivortex binding. However, 
at sufficiently large distances and thus at high tempera- 
ture it becomes negligibly small. The vortex-antivortex 
unbinding transition temperature Tkt{8) as a function of 
hole doping 8 is expected to linearly scale with the phase 
stiffness K p of the Cooper pair field. K p is linearly de- 
pendent on 8 particularly in the lightly doped region. On 
the other hand, the strength of the additional interaction 
shows a nonlinear (quadratic) dependence of 8 owing to 
the linear doping dependence of K b (~ 8) and zj(^ 8) in 
the lightly doped region 19 . 

To grasp the origin of the additional attractive in- 
teraction between vortices in a different angle, we now 
take a different procedure. Integrating over the effective 
gauge field d^ first in Eq. [2], we obtain the effective 
Lagrangian 9 of the U(l) gauge invariant particles con- 
taining the Doppler energy shift term 20 , z^Jf^d^dp, 



Z = J Di,De p e-V dT I dx2c , 
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(l/u p ,Kp,Kp), the phase 



stiffness of the Cooper pair order parameters and z^ = 
— = (z , Zj, Zj), the effective charge of the 

spinon quasiparticlcs 9 . The above Lagrangian is the low 
energy Lagrangian of the d — wave BCS theory with the 
doping dependent phase stiffness, K p (<~ 8) and the ef- 
fective charge, zj (<~ 8) 9 . 

To see the effects of the Dirac fermions on the phase 
fluctuations of the Cooper pair fields, we integrate over 



the Dirac fermion fields ignoring the local interactions of 
the Dirac fermions and the temporal fluctuations in Eq. 
[10] to find 



C = ^\V9 P - 2A\ 2 + ^;*j(V x V6 P - 2V x A) 



1 



(V x V(9 P - 2V x A). 



(11) 



The additional attractive interaction between vortices 
leads to the second term in Eq. [9] 21 as a result of the su- 
percurrent affected by the massless Dirac fermions owing 
to the Doppler shift term z^Jf^d^Op. For a brief guid- 
ance we list differences in interaction potentials between 
vortices in the s — wave and the d — wave superconduc- 
tors in Table 1 and a comparison of the d — wave BCS 
theory and the present theory in Table 2. We note that 
in the limit of q << M 2 in Eq. [7] the interaction energy 
between vortices obtained by the d — wave BCS formal- 
ism is the same as that obtained by our theory. This 
limit corresponds to the case in which the local interac- 
tions of the massless Dirac fermions are ignored in Eq. 
[10]. The additional interactions may affect the vortex 
lattice structure 22 . We believe that because our vortex 
interaction terms reveal the hole doping dependence, it 
will be of great interest in the future to examine how 
the vortex dynamics and lattice structure vary with hole 
concentration. 

We noted that the phase transition in underdoped 
cuprates falls into the XY universality class owing to 
the presence of the massive gauge field a^ which results 
from the spinon singlet pair excitations. In addition, we 
showed that the massless Dirac fermions coupled to the 
phase fluctuations of the spinon singlet pairs (leading to 
the Cooper pairs) result in the additional attractive inter- 
action of a logarithmic behavior, ^Ki,zjln\x— x'\ at short 

distances and of a power law behavior, — K ^ 2 J \x — x'\~ x 
at large distances. The attractive interaction enhances 
binding of vortex-antivortex pairs at low temperature, 
thus causing enhanced stability of the superconducting 
phase. The present study has been made based on the 
U(l) slave-boson theory concerned with the single holon 
order parameter. Thus it will be of great interest to ap- 
ply our recent SU(2) theory 23 involved with holon-pair 
boson order parameter. 
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we obtain the total interaction energy between vortices, 
\{J P v - 2d x A) li P flv (J p v - 2d x A) v , where J p y, M = 

—ii^lv^^pv ~ i'pvd^lv) is tne vortex current and P^, 
the renormalized propagator of the dual gauge field (Eq. 
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tained by the first term in Eq. [11] and the second term 
~t z J q( q +M' 2 ) ^ S e< 3 ua l to the second term in Eq. [11] in the 
limit of q « M 2 (see the text). Since we neglect the local 
interaction between Dirac fermions in Eq. [11], the second 
(logarithmic) term in Eq. [8] is not reproduced. 
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Fig. 1 The total interaction energy (solid line), logarith- 
mic interaction energy (dashed line) and additional 
interaction energy (dotted line) (in the unit of t) at 
underdoping 6 <~ 0.035 as a function of vortex dis- 
tance d (in the unit of t^ 1 ) 

Table 1 Comparison of the interaction potential between 
vortices in s — wave and d — wave superconduc- 
tors 

Table 2 Comparison of the interaction potential between 
vortices in the d — wave BCS formalism and in the 
present gauge theory formalism 



K b and K sp , are K b = ^ = 

2t Xo S and K sp = JAg = J(S - 5) 1{s . xo = J<5o' s + 205 is 
the hopping order parameter (which fits the computed re- 
sult of the U(l) slave boson mean field theory (Ref. [23])), 
where So satisfies the condition, T*(5o) = (T* : the PC 
temperature) (see the Ref. [9]). 
G. E. Volovik, Sov. Phys. JETP 58, 469 (1993). 
After integrating over the dual gauge field c M in Eq. [6], 
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TABLE I. 



s — wave superconductors 



L = 



Ldual 



-\60 v \ 



-\dxc\ 2 +icJ v (J v = dxde p ) 

V(q) = 
V(x) = K p ln\x\ 



d — wave superconductors 

^\de p \ 2 - izjj f de p + ^dipi 



-'dual 



= *k\d xc l ; 



+ icJv + jqZ 2 Jv- 



V(x) = K p ln\x\ - fz 2 j± 



- T Jv 
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TABLE II. 



d — wave BCS theory Our theory 

L = ^\89 p \ 2 - izjj f d9 p + i, n di>i L = fy\86 p - ~a\ 2 + f (d x a)-±={d x a) + \m 2 a a 2 



id^ai = 2^|9 x c\ 2 +ic.J v + ^zjJv-^Jv {Jv = dx d9 p ) L dua i = x c| 2 + icJv + i(S x c) w ^_g 2+m 2 ( a x c ) 

= £ + ^ V(q) = ^ + 

V(a;) = K p ln\x\ - f 4^ = K p ln\x\ - for q « M 2 

V(x) = (K p + ZK b zj)ln\x\ for q » M 2 
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